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We review predictions of the theory of leading twist nuclear shadowing for inclusive 
unpolarized and polarized deuteron structure functions , gf and fef and for the 
tagged deuteron structure function F^ ) (x,Q 2 ,p). We analyze the possibility to extract 
the neutron structure function F," from electron-deuteron data and demonstrate that an 
account of leading twist nuclear shadowing leads to large corrections for the extraction 
of F£ from the future deuteron collider data both in the inclusive and in the tagged 
structure function modes. We suggest several strategies to address the extraction of F^ 
and to measure at the same time the effect of nuclear shadowing via the measurement 
of the distortion of the proton spectator spectrum in the semi-inclusive eD — » e' NX 
process. We address the issue of the final state interactions in the eD — > e' NX process 
and examine how they affect the extraction of F^ ■ 

Keywords: Nuclear shadowing; Electron-Ion Collider; deuteron structure functions 
PACS Nos.:16.60.Hb, 24.85.+p, 25.30.Dh 



1. Introduction 

The nuclear shadowing phenomenon is one of key elements in the modern under- 
standing of the geometry of hadron-nucleus and nucleus-nucleus collisions. A reli- 
able theory of nuclear shadowing in hard processes is important for establishing the 
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universality of hard QCD processes. 

The initial and rather successful theory of this phenomenon was based on the 
optical (Glauber) approximation.^ However, at high energies, the picture of consec- 
utive collisions becomes inapplicable as a consequence of the increase of distances 
involved in high energy processes with increasing energiesP^El 

In addition, the eikonal (optical) approximation to inelastic processes violates 
energy-momentum conservation. In the late 60's, V. Gribov developed an approach 
to nuclear shadowing, which accounts for the specifics of high energy processes.^ 
This approach leads to the expression for the total hadron-nucleus cross section 
which is close to the Glauber approximation, but with an additional effect - inelastic 
shadowing. For the deuteron target, nuclear shadowing is unambiguously related to 
the cross section of diffraction in hadron-nucleon collisions^ 

The investigation of deep inelastic processes off nuclei is another challenge. In 
the case of scattering off the deuteron, a generalization of Gribov's ideas gave a 
possibility to evaluate nuclear shadowing in terms of diffractive parton densities of 
the nucleonP^ 

One of applications of the theory of nuclear shadowing is the possibility to 
extract the unpolarized and polarized g™ neutron structure functions from 
electron-deuteron data. 

The measurement of the non-singlet difference of the proton and neutron struc- 
ture functions, -Ff — F£ , by using the deuteron beam and spectator tagging, is one 
of the important components of the planned physics program of the electron-ion 
collider (EIC)0The main goal of the measurement is to study the flavor depen- 
dence of parton distribution functions (PDFs) in a wide kinematic region, including 
the low- a: region. In particular, the collider kinematics will allow to probe the values 
of Bjorken 10 4 for Q 2 > 1 GeV, which are a factor of ten smaller than 

in the available fixed-target data0~ln addition, the measurement of -Fjf — is a 
unique way to investigate nuclear shadowing of non-singlet (valence) PDFs in the 
non-vacuum channel. Since at small x the — F$ difference is compatible to the 
predicted nuclear shadowin g corre ction to the deuteron structure function F% (a 
few percent effect) j " ! 1 " I 1 the correct extraction of F% — F% from deuteron 
data requires an account for nuclear shadowing. We demonstrate the usefulness of 
the investigation of the deuteron tagged structure functions for the extraction of 
the nonsinglet flavor distributions. 

The accurate measurement of nuclear shadowing in nuclear PDFs and under- 
standing of its origins are of key importance for establishing the geometry of heavy- 
ion physics. This impacts the present and forthcoming RfflC data and especially 
the future LHC experiments. Therefore, it is crucial to determine the value of the 
nuclear shadowing correction, as well as its uncertainties, both to the inclusive struc- 



ture function F® and to the tagged structure function, when the spectator proton 
(a proton with momentum < 0.1 GeV/c in the deuteron rest frame) is detected 
ensuring the kinematics maximally close to the scattering off a free nucleon. 

This paper is organized as follows. In Sec. [3 we discuss our predictions for the 
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Nuclear shadowing and extraction of — FJ£ at small x from electron- deuteron collider data 3 

nuclear shadowing correction to the unpolarized F® and polarized deuteron 
structure functions. The roles of nuclear shadowing and final state interactions in 
the tagged unpolarized deuteron structure function are discussed in Sec. We 
summarize and discuss our results in Sec. 0| 

2. Leading twist shadowing and inclusive deuteron structure 
functions 

In part, this review is based on Ref. 15. For the comprehensive review of hard 
processes with the deuteron, we refer the reader to Ref. 16. The leading twist theory 
of nuclear shadowing was developed in Ref. 6 and later applied to deep inelastic 
scattering (DIS) on nuclear targets in Refs. 17-20. The theory of nuclear shadowing 
of nuclear PDFs and structure functions^ is based on Gribov's connection between 
the nuclear shadowing correction to the total hadron-deuteron cross section and 
the cross section of diffraction off a free nucleon^ Collins's factorization theorem 
for hard diffraction in dkP and QCD analyses of HERA data on hard inclusive 
diffractionP323 

On the qualitative level, Gribov's connection of nuclear shadowing to diffraction 
can be understood as a consequence of the interference between the amplitudes 
for diffractive scattering of the projectile off the proton and off the neutron of the 
deuterium target. Such interference is possible for small x, x < 5 x 10 -2 , when the 
minimum momentum transfer to the nucleon, ~ xtun, is smaller than the average 
nucleon momentum in the deuteron. 

In graphical form, the imaginary part of the forward 7*-deuteron scattering 
amplitude, which is proportional to the deuteron structure function , is presented 
in Figs, ^and [3 The dashed vertical lines denote unitary cuts placing the cut lines 
and vertices on mass-shell. Figure ^ corresponds to the impulse approximation, 
when the shadowing correction to F® is neglected. Figure|2depicts the interference 
diagram, which gives rise to the shadowing correction to F% ■ In this figure, the 
zigzag lines denote the diffractive scattering (diffractive exchange). 

One can demonstrate that the imaginary part of the interference graph de- 
creases -FrP (gives the negative shadowing correction) using the Reggeon calculus 
cutting rules of Abramovskii, Gribov and Kanchcli (AGK), which relate the shad- 
owing effects in the total and partial 

cross sections,^ S ee also Ref. 25. The top 
left graph in Fig. [21 corresponds to the diffractive final state and is proportional 
to (1 + f7 2 )(Im_4) 2 , where A is the amplitude for the photon- nucleon diffractive 
scattering and r\ = Re^4/Irru4. The top right and bottom left graphs correspond 
to the inelastic interaction with one of the deuteron nuclcons and give the com- 
mulative contribution proportional to — 4(Iul4) 2 . Finally, the bottom right graph 
corresponds to the simultaneous inelastic interactions with two nucleons and is pro- 
portional to 2(Im_4) 2 . Therefore, the sum of all four graphs in Figgis proportional 
to —(1 — 7y 2 )(Im^4) 2 = —(1 — t] 2 )/(1 + rj 2 )F^ 4 \ where we expressed the diffractive 
amplitude in terms of the nucleon diffractive structure function F®^. Therefore, 
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Fig. 1. The impulse approximation to the deuteron structure function F® ■ 




Fig. 2. All possible unitary cuts of the interference diagram giving rise to the nuclear shadowing 
correction to F® . 



the AGK cutting rules demonstrate that after summing over all final states, the 
nuclear shadowing term (the imaginary part of the interference diagram) is simply 
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proportional to the contribution of the diffractive final state (the top left graph in 
Fig.®. 

The deuteron unpolarized inclusive structure function F 2 D reada^ 

F 2 D {x,Q 2 )=F${x,Q 2 ) + F 2 n {x,Q 2 ) 

-2^—^ f ° dx P dq 2 F 2 D(4) (j3,Q 2 ,x P ,t) p D (Aq 2 t + 4{x P m N ) 2 ) , (1) 

where the first line corresponds to the impulse approximation (Fig. ^) and the 
second line corresponds to the interference term (Fig. 0). In this equation, F 2 ^ 
is the nucleon diffractive structure function; x p is the fractional loss of the proton 
longitudinal momentum; /3 ~ x/xp; xq = 0.1; q = qt + {xpmN)e z is the momentum 
transferred to the proton; pu is the deuteron charge form factor; \t\ = q 2 + (xprriN) 2 ■ 
The deuteron charge form factor can be written as an overlap of the initial and final 
state deuteron wave functions 



(2) 



p D (Aq 2 + A{xpm N ) 2 ) = / d 3 p [u{p)u{p + q) 

'3{p-(p + q)) 1 



+w(p)w(p + q) 



2p 2 (p + q) 2 2 



where u is the S'-wave component of the deuteron wave function; w is the _D-wave 
component. Note the argument of the deuteron form factor, which is a consequence 
of the correct treatment of the deuteron center of mass. Since the t-dependence 
of pu is rather moderate (compared to heavier nuclei), the integral in Eq. JQl is 
sensitive to F 2 (t) up to — t < 0.05 GeV 2 . The ratio of the real to imaginary 
parts of the diffractive scattering amplitudes rj can be calculated using dispersion 
relations over the energy or using the Regge-pole typ e p arameterization s ajF for 
the energy dependence of the diffraction cross section}^ even though the energy 
dependence differs from that for soft QCD processes, 

where ap(0) is the intercept of the effective "Pomeron" trajectory, which differs 
from th at o f the actual Pomeron, which dominates soft QCD phenomena. Using 
i] « 0.32pHlone readily observes that the correction for the real part of the diffractive 
scattering amplitude reduces nuclear shadowing by almost 20%. Note that in the 
Reggeon calculus derivation,^ it was assumed that r/ = 0, which is natural for the 
amplitudes slowly increasing with energy. This is not the case for DIS and, hence, 
the effect of rj should be taken into account. One should note that the simple final 
expression for F 2 in Eq. is due to the used closure relation for the final nuclear 
states. 

The use of the QCD factorization theorem for hard diffractioifSl allows to ex- 
tend Eq. J3> for the structure function F® to the deuteron parton distribution 
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functions"! fj/jj 

fj/D{x, Q 2 ) = fj/ p (x, Q 2 ) + fj/ n (x, Q 2 ) 

I ° dx P dq 2 fP /N ((3,Q 2 ,xjp,t) PD (4qf + 4{x P m N ) 2 ) . (4) 

Note that we use xq = 0.1 for quarks and xq = 0.03 for gluons, see the discussion 
in Ref. 18. 

The results of the calculation of the ratio of the next-to-leading order (NLO) 
structure functions / (Ff + F£) and the ratio of the NLO gluon PDFs <7£>/(2<7jv) 
are presented in Fig.|3J The solid curves correspond to Q — 2 GeV; the dashed curves 
correspond to Q = 5 GeV; the dash-dotted curves correspond to Q = 10 GeV. The 
two sets of curves for gr)/(2gjsf) correspond to the two scenarios of nuclear shadowing 
for gluons, see the discussion below and Ref. 19 for the detailed discussion. 




Fig. 3. The ratios Ff/(F^ + F 2 ™) and gn/(2gN) as functions of x. The solid curve corresponds 
to Q = 2 GeV; the dashed curve corresponds to Q = 5 GeV; the dash-dotted curve corresponds 
to Q = 10 GeV. Two sets of curves for gr)/(2g^) reflect the uncertainty in the slope of the 
t-dcpcndcnce of the gluon diffractive PDF (see details in the text). 
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As an input for our calculation, we used the HI fit to the nuclcon diffractive PDFs 
fj/N >^ where the gluon distribution was reduced by the 0.75 factor in accordance 
with the most recent HI data on hard inclusive diffraction at HERaP^ Even with 
this reduction, the gluon nucleon diffractive PDF is much larger than the quark 
diffractive PDFs. As a consequence, see Eq. (0J, we predict that shadowing in the 
gluon channel is larger than that in the quark channel and for /(F% + F^)- 

The ^-dependence of the quark diffractive PDFs fP, N was chosen to be expo- 
nential with the slopJ2SI B — 7.2 ±1.1 GeV~ 2 . Since the i-dependence of the 
gluon PDF has not been measured, the slope of the gluon PDF could be differ- 
ent. In order to take into account the uncertainty in the value of the slope B of 
the gluon diffractive PDF, we use two values: B g = 4 + 0.2 ln(10~ 3 /ir) GeV~ 2 and 
B g = 6 + 0.25 ln(10 _3 /x) GeV~ 2 . The first one corresponds to the lower end of the 
values of the J/ip photoproduction slope reported at HERA J22I while the second 
one is close to the quark slope B and to the J/ip slope reported in Ref. 30. 

The used deuteron wave function corresponds to the Paris nucleon-nucleon 
potential.^ 

Note that the results in Fig. correspond to the leading twist component of 
nuclear shadowing because they are based on the leading twist analysis of diffraction 
at HERA. When one decreases Q 2 below Q 2 > 4 GeV 2 , for instance down to Q 2 < 1 
GeV 2 , one expects a significant enhancement of the nuclear shadowing effect due 
to the enhancement of diffraction at small Q 2 by higher twist effects such as vector 
meson production. This will increase nuclear shadowing by approximately a factor 
of two P2I 

Since the diffractive structure function F% 1S known with accuracy of approx- 
imately 20% the accuracy of the calculation of the nuclear shadowing correction 
to the deuteron structure function F 2 D is 20 x 0.03 = 0.6%. Correspondingly, the 
theoretical uncertainty for the ratio of F^/F^ extracted from the deuteron data 
will be 2 x 0.6 = 1.2%, which is likely to be smaller than possible experimental 
systematic errors. 

For completeness, we also list predictions for the polarized deuteron structure 
function . Unlike the unpolarized case considered above when the shadowing 
correction was important for the extraction of F% / F£ from deuterium data because 
F% and F% are very close at small x, the polarized g\ « — at small x, which 
makes the shadowing effect a very small correction in the extraction of g" from 
polarized deuteron data, see e.g. Refs. 32 and 33. In almost complete analogy with 
the unpolarized case, the deuteron structure function can be written as a sum 
of the impulse approximation and nuclear shadowing (interference) contributions 

9?{x, Q 2 ) = (l- ^-Pd^j (<??(*, Q 2 ) + g?(x, Q 2 )) 

-2^^ p dx P dq 2 t AF D ^ (/?, Q 2 ,x P ,t) p% (4 g 2 + A(x P m N ) 2 ) , (5) 
where 1 — 3/2 Pp is the effective polarization of the proton and neutron in the 
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deuteron, which differs from unity due to the deuteron D-wave contribution (Pd = 
0.06 for the Paris nucleon-nucleon potential); p% is the electric form factor of the 
deuteron polarized in the longitudinal direction, which contains the charge and 
quadrupole form factor contributionJ^, 

AF D(4) = pDW _ F D(4) ig difference 

of the diffractive polarized nucleon structure functions. The first arrow stands for 
the helicity of the photon; the second arrow indicates the helicity of the nucleon. 
The p]j form factor has the following representation in terms of the deuteron S and 
D-w&ve components' 2 -* 1 

p% (4g 2 + 4(ipmiv) 2 ) = J d 3 p [u(p)u(p + q) 

u(p)w(p + q) { 3 (p z + q z ) 2 1\ u(p + q)w(p) ( 3 p\ 1 



.(6) 



V2 \2 (p + q) 2 2 J V2 \2p 2 2 

V(Pf(pt + qtW-(p + q)) , 3p 2 z , 3( Pz + q 



+w(p)w(p + q) 



2 p 2 (p + q) 2 Ap 2 4 {p + q) 2 

Since AF D ^ is a new and unmeasured quantity (it can be measured in polarized 
diffractive DIS on the nucleon), we cannot directly use the leading twist theory of 
nuclear shadowing to estimate the shadowing correction to gf . However, making an 
assumption that the relative strength of diffraction mediated by the non-vacuum 
exchange (responsible for the polarized structure function g\ at small x) is the 
same as that of the exchange with vacuum quantum numbers (responsible for the 
unpolarized F2 at small x), one obtains that 

AF D(4) _ p D{i) 
2 



9l t 2 



where = {g{ + g%)/2 and F 2 N = (F 2 P + F 2 n )/2. This assumption allows one 
to express the shadowing correction to the ratio of the deuteron and nucleon spin 
structure functions in terms of unpolarized diffraction on the nucleon 

9?{^Q 2 ) _ 1 _ 4 }-v 2 1 1 



2(l-|PD) 9 f l + V 2 {l-lPn)F 2 N (x,Q2) 

x [ ° dx P dq 2 t F^ {fi, Q 2 , x P ,t) p% (Aq 2 t + A(xpm N ) 2 ) . (8) 

J X 

The assumption of Eq. Q corresponds to the maximal shadowing correction. The 
additional factor of two is a source of the generic combinatoric enhancement of 
nuclear shadowing in po lari zed structure functions of few-nucleon nuclei compared 
to the unpolarized case.l^ 

Figure 01 presents the results of the calculation using Eq. (JSJ . The solid curve 
corresponds to Q = 2 GeV, the dashed curve corresponds to Q = 5 GeV and the 
dot-dashed curve corresponds to Q = 10 GeV. 

As can be seen from Fig. 0] the shadowing correction to g® could be as large as 
8% at x ~ 10 -5 . Needless to say that in order to achieve such low values of Bjorken x 
simultaneously with Q 2 > 1 GeV 2 , one needs a collider with the polarized deuteron 
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, 1.02 




Fig. 4. The ratio gf /[(l - 3/2 Po)gf ] as a function of x. The solid curve corresponds to Q = 2 
GeV; the dashed curve corresponds to Q = 5 GeV; the dash-dotted curve corresponds to Q = 10 
GeV. 



beam. The available fixed-target data® can probe only down to x ~ 0.004, 
where the shadowing correction is very small. 

It was noticed in Rcf. 16 that for a spin-one target (deuteron), the cross section 
of DIS depends on the deuteron polarization even with the unpolarized beam. The 
associated asymmetry 



T20 — 



(9) 



where a +, ° denotes the 7*-deuteron cross section and the superscript denotes the 
deuteron helicity, was estimated for x > 0.1 in the impulse approximation. Next, it 
was pointed out in Ref. 37 that nuclear shadowing in unpolarized DIS on deuterium 
leads to the values of the T20 asymmetry at the level of one percent at small- a;. 

This can be estimated as follows. The definition © allows one to immediately 
write the expression for T20 by replacing the deuteron charge form factor in Eq. 
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T2o(x, 



1 - 



F 2 D (x,Q 2 )l + V 2 



dx] 



dq^F^ (P,Q 2 : xjp,t) p 2 £ (4q 2 



4(xjpmAr) 2 )(10) 



where 



p 2 £ (Aq 2 + 4(x F m N ) 2 
u(p + q)w(p) 



= o / d "P 



^ 1 

(p- (p + q))(p- (p 



q)- 



+ w(p)w(p 4 

3pz(p z + q z )) 



-u(p)w(p + q) 



i - 



3(Pz + q z f 
(p + q) 2 

q z ? , J 



(11) 



p 2 (p + q) 2 

Note that T20 vanishes, if one ignores the D-wave component of the deuteron wave 
function or the nuclear shadowing correction. 

This effect can be also formulated in terms of the third deuteron structure func- 
tion, 6f, w hich has the following probabilistic interpretation in terms of the quark 
distribution^^ 



q°(x) + q°(x) - I {q\x) + q\x) + q-\x) + r\x)) 



(12) 



where q x is the unpolarized quark distribution in the deuteron with helicity A. The 
connection between 6f and T 20 is readily obtained using their definitions, 

\x,0 2 ^ 



b?(x,Q 2 ) 



2x 



-T 2 o(x 7 Q 2 ). 



(13) 



The factor 1 / (2x) in Eq. Ijl3(l indicates that the often discussed &f structure function 
is a rather inappropriate quantity: even small values of the physically measured T 2 q 
asymmetry correspond to huge values of fef. 

The results of the calculation of the tensor asymmetry T 2 q and the deuteron 
structure function &f are presented in Fig. |SJ The solid curve corresponds to Q = 2 
GeV; the overlapping dashed and dash-dotted curves correspond to Q = 5 GeV and 
Q = 10 GeV. 

As one can see from the left panel of Fig. the obtained T20 is at the level of 
0.6-0.7%. This agrees with the the analyses of Refs. 32 and 37. At the same time, 
6f is large at small x, as can be seen from the right panel of Fig. [3 This is a 
purely kinematic effect due to the 1/(2 a;) factor in the definition of 6f (113(1 . The 
observation of surprisingly large 6f at small x was first presented in Refs. 32, 33 
and 39. 

The HERMES mcasuremen of 6f indicates a (rapid) growth of &f when one 
decreases Bjorken x from x « 0.5 down to 1 ~ 10 -2 . However, the correspond- 
ing values of Q 2 are of the order of 1 GeV 2 and the values of Bjorken x are not 
small enough to see the predicted dramatic rise of &f towards small x. Once again, 
the study of the behavior of fof at small x will greatly benefit from the collider 
kinematics. 
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3. Nuclear shadowing and final state interactions in the tagged 
deuteron structure function 

A strategy, which is complimentary to the inclusive measurement of , is the use of 
the neutron and proton tagging. The scattering on the neutron of deuterium is then 
tagged by detecting a slow (spectator) proton. The usefulness of the tagged deuteron 
structure function for the extraction of the neutron F£ at large x was discussed in 
Ref. 41. In this work, we concentrate on the small- x region of nuclear shadowing. 
We extend the analysi^^ by taking into account the final state interactions (FSI) 
between the final nucleons. 

In the impulse approximation, the tagged deuteron structure function is given 
by the imaginary part of the left graph in Fig. [6] 

The corresponding expression is 

F 2 D (x,Q 2 ,p) =(l+i5L F?(x,Q 2 )p D (p,p), (14) 

IA \ m N J 

where p is the momentum of the spectator proton; ap (0) is the intercept of the ef- 
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Fig. 6. The impulse approximation and the final state correction to the tagged deuteron structure 
function F% . 



fective "Pomeron" trajectory; pjj(p,p) = u 2 (p) + w 2 (p), where u(p) and w(p) are the 
S'-wave and the D-wave co mpo nents of the deuteron wave function, is the unpolar- 
ized deuteron density matrbj^H (with equal momenta). The factor (l+p 2 /mAr) Qjp (°' 
comes from different invariant energies of the virtual photon-nucleus and the virtual 
photon-neutron interactions. Thus, the (1 + p z /mjv) ajp (°) factor is the flux factor 
of the interacting neutron. Note that in the derivation of the (1 + p z / m N) a,F ^ 
factor we neglected 0{p 2 /m 2 N ) and higher corrections, which are ignored in our 
non-relativistic treatment of the deuteron wave function. 

The impulse approximation receives corrections due to the final state interactions 
(the middle graphs in Fig.^J) and nuclear shadowing, which is also accompanied by 
the final state interactions (the right graph in Fig. The shaded ovals in Fig. 
denote the final state interaction. Since the FSI between two nucleons is largest at 
the small relative momentum, the FSI is accompanied by the diffractive scattering 
providing this condition. 

In the case of the tagged structure function at small spectator nucleon momenta, 
p < ^erriN, where e is the deuteron binding energy, it is legitimate to keep only the 
single and double scattering terms shown in Fig. EJ At larger spectator momenta, 
the additional contributions of triple and quadruple interactions with the target 
(not shown in Fig. IHJ are not suppressed by the p/ ^JeniN parameter and, hence, 
should be included. This will introduce a certain model dependence since those 
terms are not simply related to the elementary diffractive cross section. One should 
emphasize that this is only the case for the tagged structure function (differential 
cross section): the triple and quadruple interaction terms cancel in the inclusive 
structure function (total cross section) , which is unambiguously expressed in terms 
of the nucleon diffractive structure function, see Eq. QJ. 
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The complete expression for the tagged deuteron structure function reads 
F 2 D (x,Q\p)= + " <0> F?(x,Q 2 )p D (p,p) 

V m N J 

+ J ° dx P dq 2 F^ 4) (p,Q 2 ,x P ,t) [(2Re J d 3 p'(F(p,p') p D {p,p') 
+ ~P(p + q,p')pn{p + q,p'j) 

+ 2 J d^p'J ^WyOPi?(p'y / )V'^(p-9/2;p'-«/2))^(p-|;p' / -|)" 

- 2 1 - ^ 2 J dxpdq 2 F® {4) {(3,Q 2 ,xp,t) \F(p,p + q) p D {p,P + q) 

+ ReJ d 3 p' {F{p,p')p D {p 1 p l ) + F{p + q,p')p D {p + q 1 p l ))^ N %{p-'i-p' - |) 

+ J d 3 P 'J ^f^)^^)^ 

where iPnn is t ne continuum non-relativistic nucleon-nucleon wave function, which 
vanishes in the absence of the FSI. The arguments of V'jvw denote the final and 
initial relative momenta of the involved proton and neutron. The pd(j>,p') denotes 
the deuteron unpolarized density matrix 



The factor F(p,p') is the generalization of the nucleon flux factor discussed above 



F(p , p1 ( 1 + ^p( 1 + AV™\ (17) 

y V m N J \ m N J 

Note that the presence of the flux factor is typical for semi-inclusive cross sections. 
In total cross sections, the flux factor effects cancel in the impulse approximation 
and are of the order 0(p 2 /m 2 N ) in the interference term. Therefore, they do not 
appear in Eqs. QJ, © and l(T2")l. 

By setting ipWA = in Eq. 1)15(1 . we obtain the result of of Ref. 15. The sec- 
ond and third lines in Eq. (|15JI correspond to the FSI correction to the impulse 
approximation (the middle graph of Fig. El the rest of Eq. 1)15(1 corresponds to the 
shadowing correction, which also includes the FSI. 

In our numerical analysis of Eq. 1(15)1 . we make the following justified approxi- 
mations. First, since the dominant diffractive exchange at small x has the vacuum 
quantum numbers, the isospin of the interacting proton-neutron pair is the same as 
in the deuteron, i.e. it is zero. Note that the approximation of the exchange with 
the vacuum quantum numbers is only justified for small x, x < 10~ 2 . At larger x, 
one should also take into account transitions with the isospin change, which makes 
an estimate of the final state interactions m ore inv olved. Second, an examination of 
the isospin-zero proton-neutron phase shiftJSEl reveals that, in the appropriate 
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kinematics (small kinetic energy in the laboratory frame) , the only essential partial 
wave is the 3 Si-wave, i.e. the partial wave with L = and S — 1. Therefore, this is 
the only partial wave included in our analysis of the FSI. This allows us to suppress 
all isospin and spin indices in Eq. ((15(1 . 

The L = relative orbital momentum allows for a simple expression of the 
nucleon-nucleon continuum wave function in the momentum space 

<pNN(ki,k 2 ) = <5 3 (£i - k 2 ) + k 2 ) = S 3 (h - k 2 ) + (e2 ^°l ~ 1} , (18) 

vk i\k\\ 

where ^oi is the phase shift of the 3 5i-wave partial wave. The wave function of 
Eq. 118(1 has the usual normalization of continuum wave functions in non-relativistic 
quantum mechanics 

d% ^NNiKM i>N N (ki,k 2 ) = 5 3 (k[ - fcr) . (19) 

Since Sqi is close to 90° in the considered kinematics j ^ | 44 | ipjp^ rec j uces ^ Q the 
following simple expression, which we used in our numerical analysis 

«(fex, k 2 ) » ^tl^f 1 ^ ^ 1 ' 2 ~ |fe|2) • (20) 

Note that the continuum wave function of Eqs. (|18l) and (|20l) with <5oi ~ 90° agrees 
with the classic Bethe-Pierls result for the low relative momentum k proton-neutron 
continuum wave function in the coordinate space 

^(r)«^±^, (21) 

where cot So = —^J(.m^jk with e = 2.15 MeV being the deuteron binding energy.ESI 
Assuming that the relevant momentum k equals the root-mean-square momentum 
of the deuteron wave function, k w 130 MeV for the used Paris nucleon-nucleon 
potential, we find that the phase shift 5q ~ 110°, which is in agreement with the 
value used our analysis. 

In o ur a nalysis, we used the parameterization of <5oi provided by the SAID 
program!^! Other applications of the non-relativistic nucleon-nucleon wave function 
to estimates of the FSI can be found in Refs. 47 and 48. 

The coefficient in front of the nuclear shadowing correction term, — 2 (1 — ?7 2 )/(l+ 
t? 2 ), is dictated by the AGK cutting rules. Indeed, in the considered case of the 
proton tagging, the shadowing correction is given by the sum of the two top graphs 
in Fig.0 which enter with the coefficients 2 and —4/(1 + f? 2 ), respectively. 

Equation 1(1 5|l describes the modification of the spectrum of the produced pro- 
tons by the FSI and nuclear shadowing. We characterize the modification by the 
ratio R defined as 

Hfeq'.g- F ^ Q2 - Pl . (22) 
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The results of the numerical evaluation of the ratio R at fixed Q 2 = 4 GeV 2 and fixed 
pt and p z = as a function of x are presented in Fig. The top panel corresponds 
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Fig. 7. The ratio _R of Eq. 1151 giving the suppression of the proton spectrum by the FSI and 
nuclear shadowing. The top panel corresponds to the calculation without the FSI (tpisrN = 0); the 
lower panel is the result of the full calculation. The solid curves corresponds to pt = 0; the short 
dashed curves correspond to pt = 100 MeV; the dash-dotted curves correspond to pt = 200 MeV. 



to the calculation without the FSI [V'tvw = in Eq. (|15|) ]: the lower panel is the 
result of the full calculation. The solid curves corresponds to pt = 0; the short 
dashed curves correspond to p t — 100 MeV; the dash-dotted curves correspond to 
p t = 200 MeV. 

As one can see from the top panel of Fig. [7| nuclear shadowing (without the FSI 
effects) suppresses the spectrum of the produced protons. The suppression is larger 
for larger p t and becomes as large as 20-30% at p t = 200 MeV and x = 10~ 5 — 10~ 4 . 
This is much larger than the shadowing correction to the inclusive deuteron struc- 
ture function ■ This enhancement of nuclear shadowing is a common feature 
of semi-exclusive reactions with nuclei. Indeed, at large p t , while the impulse ap- 
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proximation term is suppressed by the nuclear wave function, the rescattering term 
becomes increasingly prominent. 

As can be seen from the bottom panel of Fig. the large effect of nuclear 
shadowing on R is mostly compensated by the final state interactions. A numerical 
analysis of Eq. 1|15[) shows that the main contribution to the FSI effect comes from 
the terms proportional to V'lvw an d that the terms containing (V']vw) 2 can be ne- 
glected. The net result is the interplay between the large and negative contributions 
of 4>nn to R, which nearly cancel each other because of the opposite signs of the 
FSI correction to the impulse approximation and the nuclear shadowing term, and 
the shadowing term without the FSI. 

From the experimental point of view, two strategies of the extraction of the 
neutron F% from the deuteron data by using the proton tagging will be possible. 
One would be to select only very low pt protons with a gross loss of statistics. Ac- 
cording to our analysis, the distortion of the proton spectrum due to the final state 
interactions and FSI effects will be minimal. The other, more promising approach 
is to measure the pt dependence of the spectrum up to pt ~ 200 MeV/c, which will 
allow to use most of the spectator protons. Provided good momentum resolution of 
the proton spectrometer, one would be able to make (longitudinal) momentum cuts 
to suppress/increase the shadowing effect and, thus, one would have an opportunity 
to independently study the enhanced nuclear shadowing. The tagged method will 
allow to extract F^ from deuteron data with the accuracy at the level of a fraction 
of percent. 

One can also use simultaneous tagging of protons and neutrons, when both 
neutron and proton are detected in the reactions — > nX and j*D — ► pX. 
In this case, the effects of nuclear shadowing and FSI will cancel in the ratio 
er 7 D ~> nX ja 1 d ^p x and the main errors in the measurement of F^ will be due 
to the determination of relative efficiencies of the proton and neutron taggers. 

One could also try to obtain the ratio Fg jF^ from the comparison of the rate 
of the tagged proton scattering events with the neutron spectator to inclusive e D 
scattering. Such a strategy could also have certain merits as it avoids the issue of 
luminosity and does not require a leading proton spectrometer. The disadvantage of 
this strategy is the sensitivity to the nuclear shadowing and FSI effects and errors 
in the acceptance of the neutron detector. One possible way to deal with the latter 
problem will be to perform measurements at very small x and large energies, where 
the ep and en cross sections are equal to better than a fraction of 1% and, hence, one 
would be able to cross-check the acceptances of the proton and neutron detectors. 

Note also that taking proton data from an independent run will potentially lead 
to another set of issues such as relative luminosity, the use of different beam energies, 
etc., which is likely to be at the level of 1%. 
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4. Conclusions and discussion 

In this brief review, we presented predictions of the theory of leading twist nuclear 
shadowing for inclusive unpolarized and polarized deuteron structure functions F® 
and and for the tensor polarization asymmetry T20. The combined role of nu- 
clear shadowing and final state interactions was analyzed for the tagged deuteron 
structure function F^(x,Q 2 ,p). The discussed effects are relevant for the collider 
kinematics of the future Electron-Ion Collider, x < 10~ 2 and Q 2 > 4 GeV 2 . 

The measurement of F® is relevant for the extraction of the neutron F^ struc- 
ture function from the deuteron data at small x. Since at small x, the proton and 
neutron F2 structure functions are close, even a small correction to F^ leads to large 
corrections to the extraction of the Ff — F^ difference for the deuteron data. We 
make predictions for the shadowing correction to F®, which enable one to extract 
FJ2/F2 from the deuteron data with an approximately 1% theoretical uncertainty. 

At moderately small x, nuclear shadowing of structure functions of light nuclei 
is dominated by the double rescattering. Thus, the theoretical analysis of F^/F®, 
where F% denotes the structure function of a nucleus heavier than deuterium, 
pre sent s another option for the extraction of the nuclear shadowing correction to 
This imposes additional constraints on the shadowing correction to F® . 

For completeness, we also consider polarized deuteron structure functions. While 
the shadowing correction to is approximately twice as large as the shadowing 
correction to F® , its effect is totally negligible in the extraction of the neutron 
<7™ from the deuteron data since g" ~ —g\ at small x and since the experimental 
errors in the measurement of g± are rather large. We find that the T20 asymmetry at 
small- x is at the one percent level so that one could try to observe it experimentally. 

The strategy complimentary to the inclusive measurement of F® is the measure- 
ment of the proton tagging. We analyzed the deuteron tagged structure function 
and showed that the combined effect of nuclear shadowing and the final state in- 
teractions only insignificantly distorts the spectrum of final protons for p t < 200 
MeV, which enables for a reliable extraction of F£ . A combined analysis of inclusive 
and semi-inclusive scattering off the deuteron coupled with a high resolution proton 
spectrometer will allow for the measurement of F% at small x with the theoretical 
uncertainty at the level of 1%. It would be a challenge to reduce the experimen- 
tal systematic errors to a comparable level. The measurement of the shape of the 
spectator spectrum would allow to determine nuclear shadowing in deuterium with 
precision by far exceeding that possible in the inclusive measurements. 

This work was supported by the German-Israel Foundation (GIF), Sofia Ko- 
valevskaya Program of the the Alexander von Humboldt Foundation (Germany) 
and the Department of Energy (USA). 
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